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ABSTRACT 
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I The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in 

^ I Stokes /, Q, and U parameters at frequencies 23, 33, 41, 61, and 94 GHz. We detect 

■ correlations between the temperature and polarization maps significant at more than 



10 standard deviations. The correlations are inconsistent with instrument noise and 



. are significantly larger than the upper limits established for potential systematic errors. 

. The correlations are present in all WMAP frequency bands with similar amplitude from 

Cf^ i 23 to 94 GHz, and are consistent with a superposition of a CMB signal with a weak 

I foreground. The fitted CMB component is robust against different data combinations 

i-^^l and fitting techniques. On small angular scales {9 < 5°), the WMAP data show the 

I ' temperature-polarization correlation expected from adiabatic perturbations in the tem- 

}^ • perature power spectrum. The data for i > 20 agree well with the signal predicted 

■ solely from the temperature power spectra, with no additional free parameters. We 

detect excess power on large angular scales {9 > 10°) compared to predictions based on 
the temperature power spectra alone. The excess power is well described by reioniza- 
^ I tion at redshift 11 < < 30 at 95% confidence, depending on the ionization history. 

' A model-independent fit to reionization optical depth yields results consistent with the 
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best-fit ACDM model, with best fit value r = 0.17 =b 0.04 at 68% confidence, including 

systematic and foreground uncertainties. This value is larger than expected given the 
detection of a Gunn-Petcrson trough in the absorption spectra of distant quasars, and 
implies that the universe has a complex ionization history: WMAP has detected the 
signal from an early epoch of reionization. 

Subject headings: cosmic microwave background, cosmology: observations, instrumen- 
tation: polarimeters 

1. INTRODUCTION 

Linear polarization of the cosmic microwave background (CMB) results from anisotropic Thom- 
son scattering of CMB photons by free electrons. By symmetry, an isotropic radiation field can 
not generate a net polarization. Any net polarization results from the quadrupole moment of the 
CMB temperature distribution seen by each scatterer. Multiple scattering suppresses polarization 
by damping the temperature anisotropy; hence, CMB polarization originates primarily from epochs 
when the opacity was of order unity or less. Standard cosmological models predict two such epochs, 
corresponding to two characteristic angular scales. The first is the decoupling surface at redshift 
z ~ 1089, when the ionization fraction abruptly falls from near unity to near zero. The acoustic 
horizon at decoupling subtends an angle ~ 1°; polarization on these scales reflects conditions in 
the photon-baryon fluid just prior to recombination. Polarization data from decoupling complement 
measurements of the temperature anisotropy. Astrophysical sources generate additional polariza- 
tion as ionizing radiation from the first collapsed objects reionizes the intergalactic medium. For 
reionization at redshift z < 50 the horizon is on large angular scales, 9 > 5°. Polarization on these 
scales directly probes the poorly-understood process of reionization. 

Since CMB polarization originates at modest opacity, the underlying temperature anisotropy 
is not heavily damped and remains observable today. Precise predictions can be made of the 
average polarization pattern expected from a given power spectrum of temperature anisotropy 
(Rees (1968); Kaiser (1983); Bond & Efstathiou (1984); Coulson et al. (1994); Kamionkowski et al. 
(1997); Zaldarriaga & Seljak (1997); Hu & White (1997); for recent reviews, see Kosowsky (1996); 
Hu (fe Dodelson (2002)). The pattern of polarization on the sky is a vector field with both an 
amplitude and direction at each point, and can be separated into two scalar fields, one giving the 
curl and the other the gradient component (called B and E modes in analogy with electromagnetic 
fields). The DASI collaboration has detected CMB polarization on angular scales ~ 0?5 (Kovac 
et al. 2002). DASI reports an E mode signal significant at 4.9cr and a TE temperature-polarization 
correlation significant at 2a. Both signals are consistent with the "concordance" cosmological model 
(spatially fiat model dominated by a cosmological constant and cold dark matter; see, e.g., Hu &i 
Dodelson (2002)) and support an adiabatic origin for the CMB temperature anisotropy. 

The Microwave Anisotropy Probe has mapped the full sky in the Stokes /, Q, and U parameters 
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on angular scales 9 > 0?2 in 5 frequency bands centered at 23, 33, 41, 61, and 94 GHz (Bennett 
et al. 2003a). WMAP was not designed solely as a polarimeter, in the sense that none of its 
detectors are sensitive only to polarization. Incident radiation in each differencing assembly (DA) 
is split by an orthomode transducer (OMT) into two orthogonal linear polarizations (Page et al. 
2003b; Jarosik et al. 2003). Each OMT is oriented so that the electric field directions accepted 
in the output rectangular waveguides lie at ±45° with respect to the yz symmetry plane of the 
satellite (see Bennett et al. (2003a) Fig. 2 for the definition of the satellite coordinate system). 
The two orthogonal polarizations from the OMT are measured by two independent radiometers. 
Each radiometer differences the signal in the accepted polarization between two positions on the 
sky (the A and B beams), separated by ~ 140°. 

The signal from the sky in each direction n can be decomposed into the Stokes parameters 

T{h) = I{n) + Q{n) cos 27 + U{n) sin 27, (1) 

where we define the angle 7 from a meridian through the Galactic poles to the projection on the sky 
of the E-plane of each output port of the OMT (Fig. 1). In principle, by tracking the orientation of 
the OMTs on the sky as the satellite scan pattern observes each sky pixel in different orientations, 
each radiometer could independently produce a map of the Stokes I, Q, and U parameters. In 
practice, the non-uniform coverage of 7 at each pixel would generate significant correlations between 
the fitted Stokes parameters, allowing leakage of the dominant temperature anisotropy into the 
much fainter polarization maps. We avoid this problem by differencing the outputs of the two 
radiometers in each differencing assembly in the time-ordered data. Denoting the two radiometers 
by subscripts 1 and 2, the instantaneous outputs are 

ATi = I{hA) + QinA)cos2^A + U{hA)sm2jA 

- I{hB) - Q(ns) cos 27B -U{hB) sin 27s (2) 

and 

AT2 = I{nA) - QifiA) cos 27A - U {ua) sin 27^ 

- I{nB ) + Q {riB ) cos 275 + U {fis ) sin 275 . 

The sum 

AT/ = ^(ATi + ATs) = I{hA) - HfiB) (3) 
is thus proportional to the unpolarized intensity, while the difference 

ATp = ^ (ATi - AT2) = Q{hA) cos 27^ + U{nA) sin 27^ - Q(nB) cos 2jb - U (ub) sin 2jb- (4) 

is proportional only to the polarization. We produce full-sky maps of the Stokes /, Q, and U 
parameters from the sum and difference time-ordered data using an iterative mapping algorithm. 
Since the polarization is faint, the Q and U maps are dominated by instrument noise and converge 
rapidly (Hinshaw et al. 2003a). 
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The Stokes Q and U components depend on a specific choice of coordinate system. For each 
pair of pixels, we define coordinate-independent quantities 

g' = Qcos(20) + C/sin(2(/.) 

C/' = C/cos(2(/))-Qsin(2(/)), (5) 

where the angle (j) rotates the coordinate system about the outward-directed normal vector to put 
the meridian along the great circle connecting the two positions on the sky (Kamionkowski et al. 
1997; Zaldarriaga k. Seljak 1997). All of our analyses use these coordinate-independent linear 
combinations of the Q and U sky maps. 

Simulations of the mapping algorithm demonstrate that WMAP can accurately recover the 
polarization pattern on the sky, even after allowing for residual calibration uncertainty in the 
individual radiometer channels. However, non-ideal instrumental signals affect the Q and U sky 
maps to a greater extent than the unpolarized I maps. The spacecraft spin about its z axis sweeps 
the beams across the sky in a direction 45° from the OMT orientation, preferentially coupling 
signals not fixed on the sky into the U map. Residual striping exists to a lesser extent in the / and 
Q maps. Systematic errors in the individual Q and U maps are not yet fully assessed; consequently, 
we defer detailed analysis of the Q or U maps to a later paper. Cross-correlations between maps are 
largely unaffected by striping or any other channel-specific signal, allowing much simpler analysis 
of the faint polarization signal than would be possible for the individual Q or U maps. This paper 
discusses the temperature-polarization (TE) correlation in the WMAP one-year sky maps. 

We compute the temperature-polarization cross-correlation using three different techniques: 
the two-point correlation function, a quadratic estimator for the power spectrum, and a "template" 
comparison in pixel space between the polarization maps and the predicted polarization given 
the observed pattern of temperature anisotropy. All three methods yield similar results despite 
disparate treatments of the data. 



2. CORRELATION FUNCTION 

The simplest measure of temperature-polarization cross-correlation is the two-point angular 
correlation function 

C^^ie) = (6) 

where i and j are pixel indices and w are the weights. To avoid possible effects of 1 // noise, we force 
the temperature map to come from a different frequency band than the polarization maps, and thus 
use the temperature map at 61 GHz (V band) for all correlations except the V-band polarization 
maps, which we correlate against the 41 GHz (Q band) temperature map. Since WMAP has a high 
signal-to-noise ratio measurement of the CMB temperature anisotropy, we use unit weight {wi = 1) 
for the temperature maps and noise weight {wj = Nj/a^) for the polarization maps, where Nj is 
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the effective number of observations in each pixel j and (Tq is the standard deviation of the white 
noise in the timc-ordcrcd data (Table 1 of Bennett et al. (2003b)). We compare the correlation 
functions to Monte Carlo simulations of a null model, which simulates the temperature anisotropy 
using the best-fit ACDM model (Spergel et al. 2003) but forces the polarization signal to zero. Each 
realization generates a CMB sky in Stokes /, Q, and U parameters, convolves this simulated sky 
with the beam pattern for each differencing assembly, then adds uncorrelated instrument noise to 
each pixel in each map. We then co-add the simulated skies in each frequency band and compute 
C^^{6) using the same software for both the WMAP data and the simulations. All analysis uses 
only pixels outside the WMAP KpO foreground emission mask (Bennett et al. 2003c), approximately 
76% of the full sky. 

Figure 2 shows C^'^{9) derived by co-adding the individual correlation functions for the fre- 
quencies 41, 61, and 94 GHz (Q, V, and W bands) least likely to be affected by Galactic foregrounds. 
The grey band shows the 68% confidence interval for the null simulations. It is clear that WMAP 
detects a temperature-polarization signal at high statistical confidence, and that signals exist on 
both large and small angular scales. We define a goodness-of-fit statistic 

= E [^MAP - iCLl)]a M„t A - {CLl)h, (7) 
ab 

where is the co-added correlation function from WMAP data, (C^-^) is the mean from the 

Monte Carlo simulations, and M is the covariance matrix between angular bins a and b derived 
from the simulations. We find = 207 for 78 degrees of freedom when comparing WMAP to 
the null model: WMAP detects temperature-polarization correlations significant at more than 10 
standard deviations. 

2.1. Systematic Error Analysis 

Having detected a significant signal in the data, we must determine whether this signal has a 
cosmological origin or results from systematic errors or foreground sources. We test the convergence 
of the mapping algorithm using end-to-end simulations, comparing maps derived from simulated 
time-ordered data to the input maps used to generate the simulated time series. The simulations 
include all major instrumental effects, including beam ellipticity, radiometer performance, and in- 
strument noise (including l/f component), and are processed using the same map-making software 
as the WMAP data (Hinshaw et al. 2003a). The Q and U maps converge rapidly, within the 30 it- 
erations required to derive the calibration solution. Correlations in the time-ordered data introduce 
an anti-correlation in the U map at angles corresponding to the beam separation, with amplitude 
0.5% of the noise in the map. This effect is independent for each radiometer and does not affect 
temperature-polarization cross-correlations. Similarly, residual l/f noise in the time series can 
create faint striping in the maps, but does not aflFect cross-correlations. 

The largest potential systematic error in the temperature-polarization cross-correlation results 
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from bandpass mismatches in the amphfication/detection chains. We cahbrate the WMAP data in 
thermodynamic temperature using the Doppler dipole from the satelhte's orbit about the Sun as a 
beam-filling calibration source (Hinshaw et al. 2003a). Astrophysical sources with a spectrum other 
than a 2.7 K blackbody are thus slightly mis-calibrated. The amplitude is dependent on the product 
of the source spectrum with the unique bandpass of each radiometer. If the bandpasses in each 
radiometer were identical, the effect would cancel for any frequency spectrum, but differences in the 
bandpasses between the two radiometers in each DA generate a non-zero residual in the difference 
signal used to generate polarization maps (Eq. 4). This signal is spatially correlated with the 
unpolarizcd foreground intensity but is independent of the orientation of the radiometers on the 
sky (polarization angle 7) . In the limit of uniform sampling of 7 this term drops out of the sky map 
solution. However, the WMAP scan pattern does not view each pixel in all orientations; unpolarized 
emission with a non-CMB spectrum can thus be aliased into polarization if the bandpasses of the 
two radiometers in each DA are not identical. This is a significant problem only at 23 GHz (K 
band), where the foregrounds are brightest and the bandpass mismatch is largest. 

We quantify the effect of bandpass mismatch using end-to-end simulations. For each time- 
ordered sample, we compute the signal in each radiometer using an unpolarized foreground model 
and the measured pass bands in each output channel (Jarosik et al. 2003). We then generate maps 
from the simulated data using the WMAP one-year sky coverage and compute C^^{9) using the 
output I, Q, and U maps from the simulation. Figure 3 shows the predicted signal at K band. 
We treat this as an angular template and compute the least-squares fit of the WMAP data to this 
bandpass template to determine the amplitude of the effect in the observed correlation functions. 
We correct the WMAP correlation functions C^^{9) and C^^ {9) at K and Ka bands by subtracting 
the best-fit template amplitudes. The fitted signal has peak amplitude of 8 jxY,? at 23 GHz and 5 
fiK? at 33 GHz. No other channel has a statistically significant detection of this effect. 

Sidelobe pickup of polarized emission from the Galactic plane can also produce spurious po- 
larization at high latitudes in the Q and U maps. We estimate this effect using the measured far- 
sidelobe response for each beam in each polarization (Barnes et al. 2003). The simplest approach 
would be to estimate the signal in each time-ordered sample, convolving the full sky sidelobe re- 
sponse with the Stokes /, Q, and U maps given the instantaneous orientation of the beams for 
each sample. Such an approach is computationally expensive. We instead approximate the signal 
in each pixel by convolving the full sky sidelobe response with the one-year Q and U maps. For 
each pixel, we fix one beam on that pixel while sweeping the other beam through all orientations 
achieved in flight. The average from the convolution yields the sidelobe contribution for the pixel in 
question. Details of the sidelobe maps are presented in Barnes et al. (2003). We correlate the side- 
lobe maps with the temperature anisotropy maps in each channel to estimate the systematic error 
in the temperature-polarization correlation. Sidelobe pickup of polarized structure in the Galactic 
plane is less than 1 /iK^ in C^^(9) at 23 GHz and below 0.1 /iK^ in all other bands. The effect of 
bandpass mismatch in the far sidelobes (as opposed to the main beam) is similarly weak, with limits 
1.3 /xK^ at 23 GHz and less than 0.05 //K^ in all other bands. We correct the polarization maps 
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for the estimated sidelobe signal and propagate the associated systematic uncertainty throughout 
our analysis. Note that all of these systematic errors depend on the Galactic foregrounds, and have 
different frequency dependence than CMB polarization. 

Other instrumental effects are negligible. We measure polarization by differencing the outputs 
of the two radiometers in each differencing assembly (Eq. 4). Calibration errors (as opposed to 
the bandpass effect discussed above) can alias temperature anisotropy into a spurious polarization 
signal. We have simulated the uncertainty in the calibration solution using both realistic gain 
drifts and drifts ten times larger than observed in flight (Hinshaw et al. 2003a). Gain drifts (either 
intrinsic or thermally-induced) contribute less than 1 /xK^ to C^^{0) in the worst band. 

Null tests provide an additional check for systematic errors. Thomson scattering of scalar 
temperature anisotropy produces a curl-free polarization pattern. A non-zero cosmological signal 
is thus expected only for the IQ (TE) correlation, whereas systematic errors or foreground sources 
can affect both the IQ and lU (TB) correlations. We also test linear combinations of radiometer 
maps which cancel the polarization signal but which test for systematic effects. We compute the 
IQ and lU correlation functions by correlating the Stokes / sum map from the Q- or V-band (as 
noted above) with the polarization difference maps [Ql - Q2)/2, {VI - F2)/2, {Wl - W2)/2, and 
(VF3 — W^4)/2. We then co-add the results with their noise weights, and compare the co-added 
result for the polarization difference maps to a similar computation for the polarization sum maps. 
The temperature (Stokes I) map in all cases is a sum map; the test is thus primarily sensitive to 
systematic errors in the polarization data. 

Table 1 shows results of the null tests. We compare C^^{9) and C^^ (9) for the sum and 
difference maps to a null hypothesis that the data consist of Stokes I and instrument noise, with no 
polarization in the Stokes Q or U maps. We break the data into 2 angular regimes to differentiate 
between signals at decoupling vs reionization. We find a clear signal detection for C^^{6) in the 
sum map for both angular scales. All other tests are consistent with instrument noise - there is no 
evidence for additional systematic errors in the temperature-polarization cross-correlation. 

2.2. Foregrounds 

Galactic emission is not a strong contaminant for CMB temperature anisotropy, but could 
be significant in polarization. WMAP measurements of unpolarized foreground emission show 
synchrotron, free-free, and thermal dust emission all sharing significant spatial structure (Bennett 
et al. 2003c). Of these components, only synchrotron emission is expected to generate significant 
polarization; other sources such as spinning dust are limited to less than 5% of the total intensity 
at 33 GHz. 

Synchrotron emission from electrons accelerated in the Galactic magnetic field is the domi- 
nant unpolarized foreground at frequencies below ~50 GHz. Although it is known to be linearly 
polarized, previous radio surveys provide little guidance for the high-latitude polarization at mm 
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wavelengths. Extrapolation of radio polarization maps (Brouw &: Spoelstra 1976) to millimeter 
wavelengths indicate a polarization fraction between 10% and 50% depending on Galactic latitude 
(Lubin & Smoot 1981). The unpolarized component has angular power spectrum q oc while 
the CMB power spectrum rises to a set of peaks on angular scales ^ ~ 1° (c/Fig 10(b) of Ben- 
nett et al. (2003c)). The angular dependence of the polarized foreground component is expected 
to be even steeper (Baccigalupi et al. 2001; Bruscoli et al. 2002; Tucci et al. 2002), suggesting 
that foreground polarization is most likely to affect temperature-polarization correlations on large 
angular scales. Radio maps at low Galactic latitude, however, demonstrate that the polarization 
intensity is not necessarily well correlated with the unpolarized intensity, complicating template 
analysis for temperature-polarization cross-correlations (Uyaniker et al. 1998, 1999). We thus use 
the frequency dependence of the measured temperature-polarization cross-correlation to separate 
cosmic from foreground signals. 

Foreground polarization above 40 GHz is faint: fitting the correlation functions at 41, 61, and 
94 GHz (Q, V, and W bands) to a single power-law C^^{e,v) = Cq^{0) {v/vqY yields spectral 
index (3 = —0.4 ± 0.4, consistent with a CMB signal (/3 = 0) and inconsistent with the spectral 
indices expected for synchrotron (/? ~ —3), spinning dust (/? ~ — 2), or thermal dust (/3 ~ 2). The 
measured signal can not be produced solely by a single foreground emission component (unless the 
fractional polarization of the foreground emission has a compensating frequency dependence, which 
seems unlikely). 

A two- component fit 



tests for the superposition of a CMB component with a single foreground component. Figure 4 
shows the resulting decomposition into CMB and foreground components. We obtain a marginal 
detection of foreground component with best-fit spectral index /3 = —3.7 it 0.8 consistent with 
synchrotron emission. We test for consistency or possible residual systematic errors by repeating 
the fit using different temperature maps and different combinations of WMAP polarization channels. 
The fitted CMB component (left panels of Fig. 4) is robust against all combinations of frequency 
channels and fitting techniques. Note the agreement in Fig. 4 between nearly independent data 
sets: the co-added QVW data (uncorrected for foreground emission) and the KKaQ data (corrected 
for foreground emission). We obtain additional confirmation by replacing the V-band temperature 
map in the cross-correlation (Eq. 6) with the "internal linear combination" temperature map 
designed to suppress foreground emission (Bennett et al. 2003c). The fitted CMB component does 
not change. We test for systematic errors by replacing the temperature map with the COBE-DMR 
map of the CMB temperature (Bennett et al. 1996), excluding any instrumental correlation between 
the temperature and polarization data. Again, the results are unchanged. 

We further constrain foreground contributions by computing the cross-correlation between 
the WMAP polarization data and temperature maps dominated by foregrounds. We replace the 
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temperature map in Eq. 6 with either the WMAP maximum-entropy foreground model (Bennett 
et al. 2003c) or a "residual" foreground map created by subtracting the internal linear combination 
CMB map from the individual WMAP temperature maps. We then correlate the foreground 
temperature map against the WMAP polarization data in each frequency band, and fit the resulting 
correlation functions to CMB and foreground components (Eq. 8). The two foreground maps 
provide nearly identical results. The fitted CMB component has nearly zero amplitude, consistent 
with the instrument noise. The fitted foreground has amplitude 0.5 it 0.1 fiK^ at vq = 41 GHz, 
with best-fit index P = —3.4 consistent with synchrotron emission. 



3. POLARIZATION CROSS-POWER SPECTRA 

In a second analysis method, we compute the angular power spectrum of the temperature- 
polarization correlations using a quadratic estimator (Appendix A). The power spectrum is the 
Legendre transform of the two-point correlation function, and is more commonly encountered for 
theoretical predictions. We compute cf^ and cf^ individually for the each WMAP frequency band, 
using uniform weight for the temperature map and noise weight for the polarization maps. We then 
combine the angular power spectra, using noise-weighted QVW data for I > 21 where foregrounds 
are insignificant, and a fit to CMB plus foregrounds using all 5 frequency bands for I < 21. Since 
foregroTind contamination is weak, we gain additional sensitivity in this analysis by using the Kp2 
sky cut retaining 85% of the sky. 

We estimate the uncertainty in each / bin using the covariance matrix M for the polarization 
cross-power spectrum. Based on our analysis of the cf^ covariance matrix (Hinshaw et al. 2003b), 
the cj^ covariance matrix has the form along the diagonal of 

M,, = <cPcf^>-<cf^>^ (9) 
^ + nTT/wi)icf'' + uee/wi) + {cf^y 

where utt and uee are the TT and EE noise bias terms, wi is the effective window function for 
the combined maps (Page et al. 2003a), cf-'^ and cf^ are the temperature and polarization angular 
power spectra, fsky = 0.85 is the fractional sky coverage for the Kp2 mask, and f^^^ = fsky/^-^^ 
for noise weighting. We take the cf"^ term from the measured temperature power spectra (Hinshaw 
et al. 2003b) and the cf^ term predicted by the best-fit ACDM model (Spergel et al. 2003) (allowing 
cf^ to vary as a function of optical depth in the likelihood analysis). Figure 5 compares the 
analytic expression for the diagonal elements of the covariance matrix to the mean derived from 
7500 Monte Carlo simulations. The analytic form (Eq. 10) accurately describes the simulations. 
We approximate the off-diagonal terms using the geometric mean of the covariance matrix terms 
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for uniform and noise weighting (Hinshaw et al. 2003b) /° 

~ ( M^M;,;, rAi. (11) 

Figure 6 shows the off-diagonal terms r^i measured from Monte Carlo simulations. The largest 
contribution, —2.8%, is at A/ = 2 from the symmetry of our sky cut and noise coverage. The 
total anticorrelation is X^A^T^o'^Ai = —0.124. Because of this anti-correlation, the error bars for 
the binned are slightly smaller than the naive estimate. A second method of estimating the 
errors relies on end-to-end simulations derived from simulated time-ordered data consisting solely of 
instrument noise (including the estimated contribution from 1/f fluctuations). We have generated 
11 "noise" sky maps each in Stokes /, Q, and U and compute the variance in TE directly from the 
variance in the simulated signal. These two approaches yield errors that are consistent to better 
than 5%. Since there are 2/ + 1 multipoles at each / value, the fractional uncertainty expected in 
the Monte Carlo variance is [2/(11(2/ + l)fsky)]^'^, in agreement with this result. 

Figure 7 shows the polarization cross-power spectra for the WMAP one-year data. The solid 
line shows the predicted signal for adiabatic CMB perturbations, based only on a fit to the measured 
temperature angular power spectrum cf^ (Spergel et al. 2003; Hinshaw et al. 2003b). Two features 
are apparent. The TE data on degree angular scales (/ > 20) are in excellent agreement with 
a priori predictions of adiabatic models (Coulson et al. 1994). Other than the specification of 
adiabatic perturbations, there are no free parameters - the solid line is not a fit to . The 
of 24.2 for 23 degrees of freedom indicates that the CMB anisotropy is dominated by adiabatic 
perturbations. On large angular scales {I < 20) the data show excess power compared to adiabatic 
models, suggesting significant reionization. 

The WMAP detection of the acoustic structure in the TE spectrum confirms several basic 
elements of the standard paradigm. The amplitudes of the peak and anti-peak are a measure of 
the thickness of the decoupling surface, while the shape confirms the assumption that the primor- 
dial fluctuations are adiabatic. Adiabatic fluctuations predict a temperature/polarization signal 
anticorrelated on large scales, with TE peaks and anti-peaks located midway between the temper- 
ature peaks Hu &; Sugiyama (1994). The existence of TE correlations on degree angular scales 
also provides evidence for super-horizon temperature fluctuations at decoupling, as expected for 
inflationary models of cosmology (Peiris et al. 2003) 

4. TEMPLATE POWER SPECTRA 

Figure 7 demonstrates that the power spectrum of temperature-polarization correlations on 
degree angular scales can be predicted using the power spectrum of the temperature fluctuations 



^'^Note that Hinshaw et al. (2003b) define off-diagonal elements in terms of the inverse covariance matrix, which 
differs from rAi by a sign. 
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alone. We use this for a third derivation of the TE cross-power spectrum, based on template 
matching in pixel space. For pixel sizes of a few degrees, the signal-to-noise ratio for the temperature 
maps is much larger than one per multipole, while the S/N ratio in the polarization maps is much 
less than one. The likelihood function for the polarization measurement then has the simple form 

log/: = (P - ^ aiPl,,^fN-\P - ^ aiP^,,^), (12) 
/ I 

where P is the measured polarization signal (a 2 A^pixei vector), ai = cf^/cf^ is the polarization 
fraction at each I, N is the pixel noise correlation matrix (a 2A'^pixei x 2iVpixei matrix) and 

Qpred{f>-) = X]"'"»(2^im("') +-2 ^^m(^^)) 
m 

Upredin) = i ^ aim{2Yim{n) --2 Yim{n)). (13) 

m 

Here ±2^im(^) are the spin harmonics, while aim are the measured coefficients for an all-sky map 
of the CMB temperature. Imposing a cut to mask the Galactic plane introduces additional corre- 
lations; we avoid this by using the "internal" linear combination temperature map (Bennett et al. 
2003c) without imposing a sky cut. 

The maps Qpred and Upred represent the predicted polarization pattern based on the observed 
pattern of temperature anisotropy. We fit these template maps to the observed Q and U polarization 
maps to derive the polarization fraction ai and thus the cf^ polarization cross-power spectrum. 
Minimizing the likelihood function yields the normal equations 

Kway = yi, (14) 

where 

yi = PN-'P^,a (15) 

and 

KW = P^redN-' Paired- (16) 

These equations show the advantages of this approach. We compare the data with a template 
in pixel space, making it straightforward to include a spatially varying noise signal. We directly 
compare the measured polarization maps to a prediction based on the measured temperature maps, 
yielding a measurement of the TE polarization cross-power spectrum in the observed sky unaffected 
by cosmic variance. We can thus more easily compute the errors on the measured polarization 
fraction. The input temperature map (Stokes /) is already corrected for foreground emission 
(much simpler in pixel space where the unpolarized foregrounds are more easily measured), greatly 
reducing the foreground contribution to the cross-power spectra. 

We thus compute the temperature-polarization cross-correlation using three disparate tech- 
niques: the two-point angular correlation function, a quadratic estimator for the power spectrum 
in Fourier space, and a template fit in pixel space. All methods are in good agreement despite their 
very dissimilar treatment of the data. All methods show a significant excess of power for I < 10. 
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5. REIONIZATION 

WMAP detects statistically significant correlations between the CMB temperature and polar- 
ization. The signal on degree angular scales {I > 20) agrees with the signal expected in adiabatic 
models based solely on the temperature power spectrum, without any additional free parameters. 
We also detect power on large angular scales (/ < 10) well in excess of the signal predicted by the 
temperature power spectrum alone. This signal can not be explained by data processing, systematic 
errors, or foreground polarization, and has a frequency spectrum consistent with a cosmological 
origin. 

The signal on large angular scales has a natural interpretation as the signature of early reion- 
ization.^^ Both the temperature and temperature-polarization power spectra can be related to the 
power spectrum of the radiation field during scattering (Zaldarriaga 1997). Thomson scattering 
damps the temperature anisotropy and regenerates a polarized signal on scales comparable to the 
horizon. The existence of polarization on scales much larger than the acoustic horizon at decoupling 
implies significant scattering at more recent epochs. 



5.1. Reionization in a ACDM Universe 

If we assume that the ACDM model is the best description of the physics of the early universe, 
we can fit the observed temperature-polarization cross-power spectrum to derive the optical depth 
T. We assume a step function for the ionization fraction Xg and use the CMBFAST code (Seljak 
Sz Zaldarriaga 1996) to predict the multipole moments as a function of optical depth. While this 
assumption is simplistic, our conclusions on optical depth are not very sensitive to details of the 
reionization history or the background cosmology. 

Figure 8 compares the polarization cross-power spectrum cf^ derived from the quadratic 
estimator to ACDM models with and without reionization. The rise in power for Z < 10 is clearly 
inconsistent with no reionization. We quantify this using a maximum-likelihood analysis 

exp(-^X^) 

|M|V2 • (17) 

Figure 9 shows the relative likelihood >C/Max(>C) for the optical depth r assuming a ACDM cos- 
mology, with all other parameters fixed at the values derived from the temperature power spectrum 
alone (Spergel ct al. 2003). The likelihood for the 5-band data corrected for foreground emission 
peaks at r = 0.17 lb 0.03 (statistical error only): WMAP detects the signal from reionization at 
high statistical confidence. 



^^Although tensor modes can also generate TE correlations at large angular scales, tensor-to-scalar ratios r large 
enough to fit the WMAP TE data are ruled out by the WMAP TT data (Spergel et al. 2003). 
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A full error analysis for r must account for systematic errors and foreground uncertainties. We 
propagate these effects by repeating the maximum likelihood analysis using different combinations 
of WMAP frequency bands and different systematic error corrections. Wc correct C^^{6) in each 
frequency band not for the best estimate of the systematic error templates, but rather the best 
estimate plus or minus one standard deviation. We then fit the mis- corrected C^^iO, z/) for a CMB 
piece plus a foreground piece (Eq. 8) and use the CMB piece in a maximum-likelihood analysis for 
r. The change in the best-fit value for r as we vary the systematic error corrections propagates the 
uncertainties in these corrections. Systematic errors have a negligible effect on the fitted optical 
depth; altering the systematic error corrections changes the best-fit values of r by less than 0.01. 

The largest non-random uncertainty is the foreground separation. We assess the uncertainty 
in the foreground separation by repeating the entire systematic error analysis (using both standard 
and altered systematic error corrections) with the foreground spectral index P = —3.7 ±0.8 shifted 
one standard deviation up or down from the best-fit value. Table 2 shows the fitted optical depth 
r and goodness-of-fit statistic for different data combinations and foreground spectral indices 
derived from the analysis of the two-point correlation function C^'^{9). The first set of rows shows 
values derived by simply co-adding the WMAP frequency channels, without any correction for 
foregrounds. Data at 41, 61, and 94 GHz (Q, V, and W bands) where foregrounds are negligible 
show similar values for r; the ~ 66 for 57 degrees of freedom indicates that the data are in 
agreement with reionized models. Adding additional low-frequency channels reduces the formal 
statistical uncertainty, but introduces non-zero foreground contamination as shown by the marked 
increase in x^- The next three sets of rows show the results when the data are separated into 
CMB and foreground components (Eq. 8). All data combinations are now in agreement; we 
obtain nearly identical values for r when fitting either the highest-frequency data set QVW or the 
lowest-frequency set KKaQ. The fitted optical depth is insensitive to the spectral index: varying the 
spectral index from -2.9 to -4.5 changes the fitted values by 0.02 or less. We adopt r = 0.17ib0.04 as 
the best estimate for the optical depth to reionization, where the error bar reflects a 68% confidence 
level interval including statistical, systematic, and foreground uncertainties. 

Spergel et al. (2003) include the TE data in a maximum-likelihood analysis combining WMAP 
data with other astronomical measurements. The resulting value, r = 0.17ib0.06, is consistent with 
the value derived from the TE data alone. The larger uncertainty reflects the effect of simultaneously 
fitting multiple parameters. The TE analysis propagates foreground uncertainties by re-evaluating 
the likelihood using different foreground spectral index. Since foreground affect only the lowest 
multipoles, the combined analysis propagates foreground uncertainty by doubling the statistical 
uncertainty in cf^ for 2 < Z < 4 to account for this effect. 

5.2. Model-Independent Estimate 

An alternative approach avoids assuming any cosmological model and uses the measured tem- 
perature angular correlation function to determine the radiation power spectrum at recombination. 
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This approach assumes that the best estimate of the three dimensional radiation power spectrum is 
the measured angular power spectrum rather than a model fit to the angular power spectrum. Given 
the observed temperature power spectrum cj"^ , we derive the predicted polarization cross-power 
spectrum cf^ (§4), which we then fit to the observed TE spectrum as a function of optical depth 
r. We obtain r = 0.16 ± 0.04, in excellent agreement with the value derived assuming a ACDM 
cosmology. We emphasize that the model-independent technique makes no assumptions about the 
cosmology. The fact that it agrees well with the best-fit model from the combined temperature and 
polarization data (Spergel et al. 2003) is an additional indication that the observed temperature- 
polarization correlations on large angular scales represent the imprint of physical conditions at 
reionization. The dependence on the underlying cosmology is small. 



5.3. Early Star Formation 

Reionization can also be expressed as a redshift Zr assuming an ionization history. We consider 
two simple cases. For instantaneous reionization with ionization fraction Xg = 1 at z < Zj., the 
measured optical depth corresponds to redshift = 17 it 3. This conflicts with measurements of 
the Gunn-Peterson absorption trough in spectra of distant quasars, which show neutral hydrogen 
present at z ~ 6 (Becker 2001; Djorgovski et al. 2001; Fan 2002). Reionization clearly did not 
occur through a single rapid phase transition. However, since absorption spectra arc sensitive to 
even small amounts of neutral hydrogen, models with partial ionization ^ 1 can have enough 
neutral column density to produce the Gunn-Peterson trough while still providing free electrons 
to scatter CMB photons and produce large-scale polarization. Direct Gunn-Peterson observations 
only imply a neutral hydrogen fraction > 1% (Fan 2002). Accordingly, we modify the simplest 
model to add a second transition: a jump from Xg = to = 0.5 at redshift Zr, followed 
by a second transition from Xe = 0.5 to Xg = 1 at redshift z = 7. Fitting this model to the 
measured optical depth yields Zr ~ 20. In reality, reionization is more complicated than simple 
step transitions. Allowing for model uncertainty, the measured optical depth is consistent with 
reionization at redshift 11 < < 30, corresponding to times 100 < U < 400 Myr after the Big 
Bang (95% confidence). 

Extrapolations of the observed ionizing flux to higher redshift lead to predicted CMB optical 

depth between 0.04 — 0.08 (Miralda-Escude 2002), lower than our best fit values. The measured 
optical depth thus implies additional sources of ionizing flux at high redshift. An early generation 
of very massive (Pop III) stars could provide the required additional heating. Tegmark (1997) 
estimate that 10~^ of all baryons should be in collapsed objects hy z = 30. If these baryons form 
massive stars, they would reionize the universe. However, photons below the hydrogen ionization 
threshold will destroy molecular hydrogen (the principal vehicle for cooling in early stars), driving 
the effective mass threshold for star formation to ~ 10^ solar masses and impeding subsequent 
star formation (Haiman et al. 1997; Gnedin & Ostrikcr 1997; Tegmark 1997). X-ray heating and 
ionization (Venkatesan et al. 2001; Oh 2001) may provide a loophole to this argument by enhancing 
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the formation of H2 molecules (Haiman et al. 2000). 

Cen (2003) provides a physically-motivated model of "double reionization" that resembles the 
two-step model above. A first generation of massive Pop III stars initially ionizes the intergalac- 
tic medium. The increased metallicity of the intergalactic medium then produces a transition to 
smaller Pop II stars, after which the reduced ionizing flux allows regeneration of a neutral hy- 
drogen fraction. The ionization fraction remains at Xe ~ 0.6 until the global star formation rate 
surpasses the recombination rate at ^ = 6, restoring Xe = 1- The predicted value r = 0.10 ± 0.03 
should be increased somewhat to reflect the higher WMAP values for the baryon density and 
normalization as (Spergel et al. 2003). The contribution from ionized helium will also serve to in- 
crease r (Venkatesan et al. 2003; Wyithe &; Loeb 2003). The WMAP determination of the optical 
depth indicates that ionization history must be more complicated than a simple instantaneous step 
function. While physically plausible models can reproduce the observed optical depth, reionization 
remains a complex process and can not be fully characterized by a single number. A more com- 
plete determination of the ionization history requires evaluation of the detailed TE and EE power 
spectra (Kaplinghat et al. 2003; Hu & Holder 2003). 



6. CONCLUSIONS 

WMAP detects statistically significant correlations between the temperature and polarization 
maps. The correlations arc inconsistent with instrument noise and are significantly larger than 
the upper limits established for potential systematic errors. The correlations arc present in all 
WMAP frequency bands with similar amplitude from 23 to 94 GHz; fitting the data to a single 
power-law in frequency yields a spectral index P = —0.4 it 0.4, consistent with a CMB signal 
(/3 = 0) and inconsistent with the measured spectral indices for Galactic foreground emission. A 
two-component fit to a superposition of CMB and Galactic foregrounds yields a positive foreground 
detection in both curl- and curl-free modes, with best-fit spectral index (3 = —3.7 ± 0.8 consistent 
with synchrotron emission of amplitude 0.5 ± 0.1 //K^ antenna temperature at 41 GHz. 

The fitted CMB component is robust against different data combinations and fitting techniques. 
On small angular scales (6 < 5°), the WMAP data show the temperature-polarization expected 
from adiabatic perturbations in the temperature power spectrum. The data for i > 20 agree 
well with the signal predicted solely from the temperature power spectra, with no additional free 
parameters. 

The data show excess power on large angular scales {9 > 10°) compared to the predictions 
based on the temperature power spectrum alone. The excess power is well described by early 
reionization at redshift Zr = 20l!^g°, corresponding to times tr = ISO^gQ^ Myr after the Big Bang 
(95% confidence). A model-independent fit to reionization optical depth yields results consistent 
with the ACDM model. Our best estimate for the optical depth isr = 0.17ib0.04 (68% confidence) 
where the error terms include statistical, systematic, and foreground uncertainties. This value is 
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larger than expected given the detection of a Gunn-Peterson trough in the absorption spectra of 
distant quasars, and imphes that the universe has a complex ionization history. 

The WMAP detection of early reionization opens a new frontier to explore the universe at 
redshift 6 < z < 30. WMAP''s sensitivity to reionization is currently limited by instrument noise, 
both as direct statistical uncertainty and in the ability to better model and remove faint polarized 
foregrounds. Instrumental effects do not limit analysis of temperature-polarization correlations. 
The TE power spectrum and covariance matrix are available at http://laiiibda.gsfc.nasa.gov. 
We are currently performing a more complete set of systematic error analyses in the individual Q 
and U maps. A future data release will include full-sky polarization maps and polarization power 
spectra. 

The WMAP mission is made possible by the support of the Office of Space Sciences at NASA 
Headquarters and by the hard and capable work of scores of scientists, engineers, technicians, 
machinists, data analysts, budget analysts, managers, administrative staff, and reviewers. 
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A. Quadratic Estimator for Temperature-Polarization Power Spectrum 

We estimate the temperature-polarization power spectrum from pixelized sky maps using the 
following formalism. We begin by expanding the temperature and polarization fluctuations in 
generalized spherical harmonics 

T{n) = J2(^lmYim{n) (Al) 

Im 

Q{h) ± iU (ra) = a^2,im T2^m(n) (A2) 

Im 

We then decompose the polarization fluctuations into E and B like pieces 

0'±2,im = Elm =•= iBi^. (A3) 
We can use the basic properties of the spherical harmonics 



I' 



dn NYUfi)NYirm,{n) = ,5fc' (A5) 



to derive 



Eim = \j dn [ Q{n) ( 2YCm{n) + -a^^C")) 
-iU{n){2YCm{h)- -2^4 (n)) ] 
---\jdn [ [/(n) ( ^Y^h) + -2Yah)) 

+i Q{n) ( 2YCr^{h) - -2YCm{n)) ] • (A6) 



We can now generalize the approach of Hivon et al. (2002) to estimate the coupling terms. We 
multiply the temperature and polarization maps by a weighting function 



Tirn = J dnw' (n) T{h)Yil{h) (A7) 

Elm = IJ dnw^'ifi) [Q{n){2YiUn)+ -2YCm{n)) 

-iU{n) ( 2YL{n) - -2YC^{n))] (A8) 

Bim = -\j dnwP{n) [U {h) { 2YC^[n) + -2l^4(^)) 

+iQ{h) ( 2Yi*m{n) - -2YCm{n))] ■ (A9) 
We expand the weighting function in spherical harmonics 

"^(^) = 'Y'^lmYlra{n), (AlO) 
Im 
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and combine with equations Al - A3 to yield 



Tim 
Elm 



XI '^J"m"Tl'm' I dhYiira'{n)Yitira''{n)Y{^{h) 
I'm'V'm" 



in) ( 2Y,m'{n) 2YCm{n) + -2ilw(n) -2YCm{n)) 



Im 



I'm'V'm" 

+iBlm J dh Yif'm"{n)i 2Yi'm'in) 2>"4("') " -2>l'm'(^) -2Yi*^ih)) 

\ X '^^'m" ^I'm' [ dhYi>rm''in)i2Yi'rn'in) 2Yi*^{n)+ -2Yi'm'{n) -2Yi*^{n)) 

I'm'V'm" 

-lEim J dh Yi»^„{h) ( 2Yi'm'{n) 2Yi^{h) - -2Yi'm'in) -2^^(»^)) 



(All) 



We can then use 

d'Tl AT^ml"-) N'Yi'm'ih) N"Yv'm"in) 



I' 



\N+m 



(2Z + l)(2r + l)(2Z" + l) 



1/2 



to compute 



7.TB 
^EE 



I r I" 

-N N' N" 



( cr \ 

JTE 
H' 

H' 



I I' I" 



-m m' m" 



(A12) 



v 



EE 



(A13) 



These expressions can be reduced using the symmetry and orthogonahty properties of S-j symbols, 
as given in Eqs. 1.8 and 1.14 of Rotenberg et al. (1959). In particular, imaginary terms drop out, 
and summations over products of 3-j symbols with — m, m/ and m" in the bottom row evaluate to 
1/(2/" + 1). After some algebra, the coupling terms reduce to 



j^TT,TT 



^TE,TE 



][/fEE,EE 



-^^1.^1" 



j^TB,TB 



87r 



v 



pp 



v 







I I' I" 

-2 2 



/ V I" \ I I' I" 
-2 2 J \ 2 -2 



(A14) 



(A15) 



+ 



I I' I" 
2-2 
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where 



/ V I" 
-2 2 



+ 



j^BB,EE 



{21' + 1 




I 


v 


I" 


2 


-2 





I 


/' 


r' 


-2 


2 





I 


Z' 


1" 


2 


-2 






m 



(A16) 




(A17) 



(A18) 



with a and 6 referring to either T or P. All of the other coupling terms are zero. Note that if 
we use different weighting functions for T, Q and U, we increase the coupling between E and B 
modes. 



B. Uniform Temperature Weighting 

If we use the full sky to compute the temperature spherical harmonic terms, then the cross- 
correlation term and its error matrix becomes particularly simple. For this case, Wqq = I/a/Itt and 
all other coupling terms are 0. In this limit, the measured cj!^ is just a constant times the true cf^ 

= -f (Bl) 

where 

/ = / -M^)^ (B2) 
The covariance matrix for these terms are diagonal. 
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Table 1: Null Tests for IQ and lU Sum and Difference Data^ 

Sum Map^ Difference Map^ 



Correlation Range DOF Probability'^ Probability^ 



IQ 


9 <5° 


20 


62.1 


3x10"^ 


23.6 


0.26 


IQ 


e>h° 


58 


145.1 


2x10"^ 


66.0 


0.22 


lU 


e<5° 


20 


30.9 


0.06 


10.8 


0.95 


lU 


e>b° 


58 


66.1 


0.22 


50.4 


0.95 



"x^ comparison of the WMAP correlation functions C'^{9) and C'^ {6) to a null hypothesis of CMB temperature 
anisotropy and instrument noise, but no polarization. Temperature-polarization signals of cosmic origin should 
contribute to only C"^{9) in the sum maps. All other tests are consistent with the null hypothesis. 

'Polarization sum maps (Ql+Q2)/2, (Vl+V2)/2, (Wl+W2)/2, and (W3+W4)/2 co-added with noise weights. 

^Polarization difference maps (Ql-Q2)/2, (Vl-V2)/2, (Wl-W2)/2, and (W3-W4)/2 co-added with noise weights. 

"^Probability to randomly obtain larger than the measured value. 
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Table 2: Rcionization Optical Depth*^ 
Data Set Method P ^ f{>WMAPy 



vw 


Co-Add 


^•^^-0.03 


67.0 


0.159 


QVW 


Co-Add 


0.15 ± 0.04 


66.2 


0.176 


KaQVW 


Co-Add 


0.14 ±0.03 


97.1 


0.001 


KKaQVW 


Co-Add 


0.30 ± 0.02 


359.8 


0.0 


KKaQ 


Co-Add 


0.29 ± 0.01 


476.6 


0.0 



QVW 


Fit 


-2.9 


12+^-'^^ 

^■^■^-0.08 


65.2 


0.201 


KaQVW 


Fit 


-2.9 


20+°-^^ 


69.8 


0.101 


KKaQVW 


Fit 


-2.9 


0.22 ± 0.04 


60.9 


0.313 


KKaQ 


Fit 


-2.9 


0.20 ± 0.04 


58.7 


0.404 



QVW 


Fit 


-3, 


.7 


1 

"^-0.07 


66. 


.1 


0. 


,180 


KaQVW 


Fit 


-3, 


.7 


19+013 

'-'•-'-^-0.06 


68. 


.9 


0. 


,117 


KKaQVW 


Fit 


-3, 


.7 


0.17 ±0.03 


55. 


.4 


0. 


,491 


KKaQ 


Fit 


-3, 


.7 


0.18 ± 0.04 


48. 


.0 


0. 


,772 


QVW 


Fit 


-4, 


,5 




66, 


,6 


0. 


,170 


KaQVW 


Fit 


-4, 


.5 


1 5+°-l'l 
'-'•-'-'3-0.05 


68. 


.2 


0. 


,140 


KKaQVW 


Fit 


-4, 


.5 


0.16 ±0.03 


57. 


.8 


0. 


,419 


KKaQ 


Fit 


-L 


,5 


O.IG ± 0.01 


51, 


,1 


0. 


,G51 



"Optical depth r fitted from C^^(6>) for various combinations of data and foreground corrections in a ACDM cosmology. 

There are 57 degrees of freedom for each fit. 
''68% confidence statistical uncertainties 

'^Fraction of 1000 simulations of reionized ACDM models with larger than WMAP value. 
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Fig. 1. — Geometry for Stokes Q and U parameters. WMAP measures polarization by differencing 
two orthogonal polarization channels, then solving for Q and U as the spacecraft compound spin 
projects the OMT onto the sky at different angles 7 relative to the Galactic meridians. All analysis 
uses coordinate-independent quantities Q' and U' defined with respect to the great circle connecting 
a pair of pixels (see text). 
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Fig. 2. — Temperature-polarization correlation function for WMAP co-added QVW data. The 
gray band shows the 68% confidence interval for similar co-added data taken from Monte Carlo 
simulations without polarization. The inset shows data for 9 < 10°. The data are inconsistent with 
no temperature-polarization cross-correlations at more than 10 standard deviations. Note that the 
data are not independent between angular bins. 



-27- 




D 30 60 90 120 150 1# 
4ligplar Separation (Degrees) 




Fig. 3. — Angular templates for potential systematic errors caused by bandpass mismatch between 
the 2 radiometers in each differencing assembly. We fit this template to the correlation functions 
from each DA to detect or limit systematic errors related to bandpass mismatch in the main beam. 
The effect is significant only in K and Ka bands, which have the brightest unpolarized foregrounds. 
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Fig. 4. — Fitted CMB (left) and foreground (right) components from a multi-frequency decompo- 
sition of the measured two-point correlation functions. Top panels show the IQ (TE) correlation, 
while bottom panels show lU (TB). The CMB component is shown in units of thermodynamic tem- 
perature, while the foreground is shown in antenna temperature evaluated at 41 GHz. Different 
colors show the effect of using different temperature maps in the cross-correlation, or including dif- 
ferent polarization frequency channels in the CMB-foreground decomposition. "Co-Add" refers to 
a noise- weighted linear combination of the correlation functions computed for individual frequency 
channels. "Fit" refers to a 2-component fit (Eq. 8) using the specified polarization frequency chan- 
nels. The grey band shows the 68% confidence interval for the CMB component for the KKaQVW 
fit (which has the smallest statistical uncertainty) assuming CMB temperature anisotropy and 
instrument noise, but no CMB polarization. "Combination" and "COBE-DMR" replace the tem- 
perature map in Eq. 6 with maps with reduced foreground emission: either the WMAP internal 
linear combination map or the COBE-DMR map of the CMB temperature. "MEM Model" and 
"ILC Residual" replace the temperature map in Eq. 6 with maps dominated by foreground emis- 
sion: either the WMAP maximum-entropy foreground model or the residual map produced by 
subtracting the internal linear combination map from the individual temperature maps at each 
frequency. The fitted CMB component is stable as different frequency channels and data sets are 
analyzed. Foreground emission is faint compared to the cosmic signal. 
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Fig. 5. — Diagonal elements of the covariance matrix for the cf^ polarization cross-power spectrum. 
Points show the diagonal elements computed from 7500 Monte Carlo simulations. The solid line 
shows the analytical model (Eq. 10). Note we multiply M// by (^) to match the units in Figures 
7 and 8. 
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Fig. 6. — Off-diagonal correlations vai in the covariance matrix for the cf^ polarization cross- 
power spectrum, computed from simulations. All values are normalized to r^i = 1 at AZ = 0. The 
dotted line shows r^i = for comparison. The anti-correlation at AZ = 2 results from the spatial 
symmetry of the sky cut and noise coverage. 



-31 - 



-3. 



o 



+ 



4 



2 





■1 



Reionization 



Adiabatic Prediction 








100 



200 300 

Multipole I 



400 



500 



Fig. 7. — Polarization cross-power spectra c^^ for the WMAP one-year data. Note that we plot 
(Z -I- l)/27r and not l{l + 1)/2'!t cf^. This choice emphasizes the oscillatory nature of cj^. For 
clarity, the dotted line shows q = 0. The solid line is the predicted signal based on the cj'^ power 
spectrum of temperature anisotropy - there are no free parameters. The TE correlation on degree 
angular scales (/ > 20) is in excellent agreement with the signal expected from adiabatic CMB 
perturbations. The excess power at low I indicates significant reionization at large angular scales. 
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Fig. 8. — WMAP Polarization cross-power spectra cj^ (filled circles) compared to ACDM models 
with and without reionization. The rise in power for Z < 10 is consistent with reionization optical 
depth r = 0.17 lb 0.04. The error bars on WMAP data reflect measurement errors only; adjacent 
points are slightly anti-correlated. The grey band shows the 68% confidence interval from cosmic 
variance. The value at Z = 7 is particularly sensitive to the foreground correction. 
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Fig. 9. — Likelihood function for optical depth r for a ACDM cosmology, using all 5 WMAP 
frequency bands fitted to CMB plus foregrounds with foreground spectral index /? = —3.7. After 
including systematic and foreground uncertainties the optical depth is consistent with a value 
T = 0.17 with 95% confidence range 0.09 < r < 0.28. 



